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Abstract 

A simple model wavefunction, consisting of a linear combination of two free-particle Gaussians, 
describes many of the observed features seen in the interactions of two isolated Bose-Einstein 
condensates as they expand, overlap, and interfere. We show that a simple extension of this idea 
can be used to predict the qualitative time-development of a single expanding BEC condensate 
produced near an infinite wall boundary, giving similar interference phenomena. We also briefly 
discuss other possible time-dependent behaviors of single BEC condensates in restricted geometries, 
such as wave packet revivals. 
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I. INTRODUCTION 



It can be argued that much of the early success of quantum theory can be traced to the 
fact that many exactly soluble quantum models are surprisingly coincident with naturally 
occurring physical systems, such as the hydrogen atom and the rotational/ vibrational states 
of molecules. Many other exemplary quantum mechanical models, which have historically 
been considered as only textbook idealizations, have also recently found experimental real- 
izations. Advances in areas such as materials science or laser trapping and cooling of atoms 
have allowed the production of approximations to a number of systems which had typically 
been relegated to lists of pedagogical examples. While some such examples have found use 
as devices, many others have been applied to the study of fundamental quantum behavior. 

For example, one-dimensional quantum wells (infinite and finite) are a staple of textbooks 
and have found use in modeling quantum dots and other structures [p.], including the use 
of asymmetric wells [2]. Two-dimensional quantum mechanical 'standing waves' have been 
observed in a wide variety of geometries p|, while evidence for bound quantum states of 
the neutron in the Earth's gravitational field 0] (a problem which is often described as 
the 'quantum bouncer' in the pedagogical literature) has recently been presented. The 
"generation of nonclassical motional states", such as coherent and squeezed states [^J of 
harmonically trapped ions, has also been demonstrated and studied in detail as an example 
of time-dependent designer wave functions in the most familiar of all model potentials. 

The experimental realization of Bose-Einstein condensates (BECs) |6| has allowed for an 
even wider variety of fundamental tests, including the " Observation of interference between 

Q Q 

two Bose condensates" y\. In the original experiment [7|, two samples of sodium atoms were 
evaporatively cooled ll well below the transition temperature to obtain condensates 11 such that 
they were initially well-separated in a double-well potential. The two condensates were 
then allowed to freely expand and interference effects were observed in the resulting overlap 
region, while no similar effects were seen for a single expanding condensate. Similar effects 
have been seen in other experimental realizations [S| and more recently have been observed 
with up to 30 uncorrelated condensates [l| produced in an optical trap. 

As we will briefly review in Se c. ITT1 a simple model wavefunction consisting of a linear 
combination of two [l0|] (or more |9|) Gaussian terms, ipia){x : t; Xq), (one for each conden- 
sate) captures many of the salient features observed experimentally. A (single particle) 
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wavefunction of the form 



i[>(x, t) = N \_ip(G) (x, t- x A ) + e^V(G) (x, t; x B )] (1) 



can be used, with x = xa .xr describing the locations of the two isolated condensates and 
with fixed relative phase (and a normalization constant N) and we will discuss the 



predictions of this simple model in the next section. 

Such linear combination solutions (especially of Gaussians) have been frequently used in 
the pedagogical literature to describe the time-development of wave packet solutions of the 
ID Schrodinger equation (SE) describing an otherwise free-particle impinging on an infinite 
wall or barrier (or "free particle on the half-line".) "Mirror" or "image" solutions of the 
form 

j N[i/;(x,t)-i>(-x,t)} forx<0 
ip(x,t) = < (2) 
I for < x 

(for a particle restricted to x < 0) satisfy the ID free-particle SE (if ip(x,t) does) and also 
automatically satisfy the appropriate boundary condition at the infinite wall (assumed to 
be at x = 0) for all times. Such solutions have been used in a variety of pedagogical applica- 
tions Q| - 3], but also in a research context to discuss the deflection of ultracold quantum 
particles (wavepackets) from impenetrable boundaries or mirrors Q]. Such analyses natu- 
rally explain the spatially oscillatory behavior of the position-space probability density as 
the wave packet 'hits' the wall (as observed in numerical calculations) as the interference 
between two overlapping terms, much like the observed BEC effect. 

In this Letter, we note that a single BEC condensate, produced near an infinite wall 
boundary and allowed to expand freely, will likely exhibit interference effects describable 
by such 'mirror' or 'image' solutions as in Eqn. (J2J), and we discuss this in Sec. IIHI We 
also extend such ideas to a localized BEC produced between two infinite boundaries (as in 
an infinite well potential) to very briefly discuss other possible effects, such as wave packet 
revivals. 



II. SIMPLE MODEL OF INTERFERING BOSE-EINSTEIN CONDENSATES 

A description of the initial (single particle) one-dimensional wavefunction for two sep- 
arated BEC condensates using isolated Gaussian forms has been made by the authors of 
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Ref. lOj . who then analyze the resulting time-development of the BECs using a Wigner 
distribution approach. (A generalization to multiple BECs is given in Ref. As a review 
of such an approach and its successes in qualitatively modeling more sophisticated analyses 
(and the experimental data), consider two BEC condensates, separated by a distance d and 
initially centered at Xq = ±d/2, described by the time-dependent free-particle wavefunction 

1P2BEC(X, t) = = - \ e -i*-m 2 /W\Wt/t ) + e#e -(x+d/2)V2^(l+^o)l (3) 

vVW + it /to) 1 J 
where the normalization factor is given by 

N = JL (l + cos(4>)e- d2/ ^y 1/2 . (4) 



N 2 

P2BEC[X,t) 



(6) 



V2 

The time-dependent spatial width for a single such Gaussian term is given by 

As t = A = ^/i + ( t/to )2 where t = ^f (5) 

and the corresponding spread in momentum-space (again, for a single Gaussian) is Ap t = 
Ap = h/{(3\/2). The time-dependent probability density for the two condensate state is 
then given by 

vnpt L V toPt 

where the cross-term describes the interference effect. 

For each Gaussian contribution, there are momentum components of order p ~ Ap Q = 

h/ (3\f2 so that the time (To) it takes such components to drift from one condensate and 

overlap with the other is of order To ~ d/{jp/m) ~ \/2dmf3/h. For condensates which 

are initially highly localized and well-separated, the time dependent position-spread is then 

dominated by the (t/to) 2 term since (To/to) 2 ~ 2(d/ /5) 2 >> 1. In that limit, the oscillatory 

term is then approximately given by 

/ xdt\ ( xdt \ 

cos 10 + -yr- — > COS 10 + —zrr (7) 



P?t j v P 2 t 

so that the local wavelength variations seen in the interference pattern are time- dependent 
and scale like 

2tt xdt ht 

T x = kx = W or A = ^' (8) 

just as in Eqn. (1) of Ref. [7j. The time-dependent real and imaginary parts of the in- 
dividual components of this simple wavefunction are nicely consistent with more detailed 



calculations [8| where the BEC is 11 characterized by a phase that varies quadratically ...across 
the condensate 11 ; for example, compare the pedagogical illustration of an expanding p = 
Gaussian wavepacket in Fig. 2 of Ref. [16] with Fig. 1 of Ref. jjj- Finally, the corresponding 
momentum-space probability density is given by 

_ 4AT^ fpd\ v 



\<p2BEc( P ,t)\ 2 = ( ^ ) e-*r (9) 



where a = /3/h so that there is indeed structure in momentum space at integral multiples of 
p — h/d, as discussed in more detail in Ref. Q|. Thus, in many important ways, the simple 
wavefunction in Eqn. (J3J) encodes much of the physics observed in the interference of two 
expanding BEC condensates. 

III. SINGLE BEC NEAR A INFINITE WALL BOUNDARY AND RELATED EF- 
FECTS 



Q, 



Motivated then by earlier pedagogical papers on 'mirror' or 'image' solutions 
we can imagine a single BEC condensate produced close to an infinite barrier. If the barrier 
is located at x = and the single condensate is produced at x = —d/2, the resulting 'mirror' 
solution in Eqn. (j2J), for x < 0, can be described by the form in Eqn. Q with <p = tt (so 
that cos(0) = —1) and with normalizations simply related by N = y/2N. The resulting 
time-dependent solution will then exhibit the same type of interference patterns observed 
for two isolated condensates. One possibility for such an infinite barrier might be an atomic 
mirror of the type successfully used in a number of atomic physics applications 

The addition of a second infinite wall barrier (say at x — —d) to such a case might then 
be modeled by the standard infinite well problem of textbooks. The time development of 
wave packet propa gati on in this system can then be described in terms of an infinite number 
of image solutions |20| and discussions related to this approach go back to at least Einstein 



and Born j21|. For a single BEC condensate in this restricted geometry, modeled as a p = 
Gaussian, in addition to the spreading/coherence time (to) and the time to overlap the other 
real or image condensate (To), the only other relevant time scale is the quantum revival time 



22j T rev . For the quantized energy eigenvalues in an infinite well of width d (as imagined 
here) the revival time for an arbitrary localized wave packet is T rev = 4md 2 /h7r j^]. The 
ratio of revival time to overlap time is then T rev /To ~ 2djnI\xQ >> 1. In the original two 
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BEC experiment |]J, the two condensates are allowed to fall freely as they expand, so a 
more detailed analysis of a specific experiment realization would be required to determine if 
the revival time is too long to be observed. One should note, however, that for the special 
case of a p = Gaussian waveform produced precisely in the center of such an infinite well, 
because only even energy eigenstates are excited, the effective revival time is actually T rev /8 



22j for this very special geometry, which gives almost an order-of-magnitude shorter fall 
time to achieve a revival. 

Other arrangements of two infinite plane barriers can also be imagined to give rise to 
interesting BEC interference effects, which can be modeled using 'mirror' or 'image' methods. 
For example, two such infinite walls can be placed at right angles, to form a 'corner (90°) 
reflector', defined by the potential 

{0 for < x and < y 
. (10) 
+oo otherwise 

A solution of the form 

VWnerO, y,t) = N [i/j(x, y; t) - if)(-x, y, t) - ip(x, -y; t) + 4>(-x, -y; t)} (11) 

making use of three auxiliary 'image' components, solves the Schrodinger equation in the 
allowed region for any free-particle solution if)(x,y;t), as well as satisfying the boundary 
conditions at the two walls. The normalization factor can be obtained explicitly in the case 
of an Gaussian solution, as above. The same construction can also be employed for other 
angles between the two walls, O, using familiar examples from optics, such as for the cases 
of = 45° and 60°. 
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